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The formation mechanism of brown dwarfs 
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ABSTRACT 

We present results from the first hydrodynamical star formation calculation to demon- 
strate that brown dwarfs are a natural and frequent product of the collapse and 
fragmentation of a turbulent molecular cloud. The brown dwarfs form via the frag- 
mentation of dense molecular gas in unstable multiple systems and are ejected from 
the dense gas before they have been able to accrete to stellar masses. Thus, they can 
be viewed as 'failed stars'. Approximately three quarters of the brown dwarfs form 
in gravitationally-unstable circumstellar discs while the remainder form in collapsing 
filaments of molecular gas. These formation mechanisms are very efficient, producing 
roughly the same number of brown dwarfs as stars, in agreement with recent observa- 
tions. However, because close dynamical interactions are involved in their formation, 
we find a very low frequency of binary brown dwarf systems 5%) and that those bi- 
nary brown dwarf systems that do exist must be close ^ 10 AU. Similarly, we find that 
young brown dwarfs with large circumstellar discs (radii <; 10 AU) are rare (« 5%). 

Key words: accretion, accretion discs - binaries: general - brown dwarfs - circum- 
stellar matter - hydrodynamics - stars: formation - stars: mass function. 



1 INTRODUCTION 

The existence of brown dwarfs was incontrovertibly demon- 
s trated for the first t ime by the discovery of Gliese 229B 
(Nakajima et al. 1995), a cool brown dwarf orbiting an M- 
dwarf. In the same year, other candidates later confirmed 
to be free-floating brown dwarfs were announced (e.g. Teide 
1 by Rebolo, Zapatero-Osorio & Martin 1995), along with 
PP1 15 which was later discovered to be a binary brown 
dwarf (Basri & Martin 1999). Observations now suggest that 
brown dwarfs are as common as stars, although stars domi- 
nate in terms of mass (e.g. Reid et al. 1999). 

Despite the abundance of brown dwarfs, their formation 
mechanism is currently a mystery. The typical thermal Jeans 
mass in molecular cloud cores is » 1 Mq (Larson 1999 and 
references therein). Thus, the gravitational collapse of these 
cores might be expected to form stars, not brown dwarfs. 

There are two obvious routes by which brown dwarf 
systems (i.e. brown dwarfs without stellar companions) may 
form. First, they may result from the collapse of low-mass 
cores (masses ^0.1 Mq) that are smaller (radii 0.05 pc) 
and denser (n(H_2)<; 10 7 cm -3 ) than the cores that are typ- 
ically observed (i.e. they have low masses yet are still Jeans 
unstable). Thus, brown dwarfs would be 'low-mass stars'. 
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Such low-mass bound cores have not yet been observed, but 
they would be difficult to detect because of their small sizes 
and low masses, and rarely detected due to their short dy- 
namical time-scales (~ 10 4 years). Observations are begin- 
ning to reach this mass regime (e.g. Motte, Andre & Neri 
1998), although the low-mass clu mps found thus far prob a- 
bly are not gravitationally bound (Johnstone et al. 2000). 

The second possibility is that brown dwarfs form in 
higher-mass cores but are prevented from accreting enough 
mass to exceed the hydrogen-burning limit. If such a core 
fragments to form an unstable multiple system, this may 
be achieved by the dynamical ejection of a fragment from 
the core, cutting it off from the reservoir of gas, and thus 
preventing it from accreting to a stellar mass. In this case, 
brown dwarfs would be 'failed stars'. This ejecti on m echa- 
nism has been propos ed by Reipurth & Clarke (2001) and 
Watkins et al. (Il998b|). 



In this letter, we present results from the first hydro- 
dynamical calculation to demonstrate that a large number 
of brown dwarfs can be formed during the fragmentation 
of a molecular cloud. All the brown dwarfs are formed by 
the ejection of fragments from unstable multiple systems. 
In section 2, we briefly describe the numerical method and 
the initial conditions for our calculation. In section 3, we 
present results from our calculation and compare them with 
observations. Finally, in section 4, we give our conclusions. 
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2 COMPUTATIONAL METHOD AND INITIAL 
CONDITIONS 

The calculation presented here was performed using a three- 
dimensional, smoothed particle hydrodynamics (SPH) code 
based on a version originally developed by Benz (Benz 1990; 
Benz et al. 1990) . The smoothing lengths of particles vary in 
time and space, such that the number of neighbours for each 
particle remains approximately constant at iV nc igh = 50. We 
use the standard form of artificial viscosity (Monaghan & 
Gingold 1983) with strength parameters a v = 1 and Pv_~ 2. 



Further details can be found in Bate, Bonnell & Price ( 1995 ) 
The code has been parallelised by M. Bate using OpenMP 



2.1 Opacity limit for fragmentation and the 
equation of state 

When the collapse of a molecular cloud core begins, the 
gravitational potential energy released easily radiates away 
so that the collapsing gas is approximately isothermal (e.g. 
Larson 1969). Thus, the thermal pressure varies with den- 
sity p as p cx p v where the effective polytropic exponent, 
r\ = dlog[p]/dlog[p] ~ 1. Fragmentation is allowed because 
the Jeans mass decreases with increasing density if r\ < 4/3. 

The opacity limit for fragmentation occurs when the 
rate at which energy is released by the collapse exceeds 
the rate at which the gas can cool (Low & Lynden-Bell 
1976; Rees 1976). The gas then heats up with rj > 4/3, the 
Jeans mass increases, and a Jeans-unstable collapsing clump 
quickly becomes Jeans-stable so that a pressure-supported 
fragment forms. The gas b egins to heat significantly a t a 



density of « 10 3 g cm 3 (Masunaga fc Inutsuka 200C) 



The pressure-supported fragment initially conta ins sev- 
eral J upiter-masses (Mj) and has a radius of « 5 AU (Larson 



196S). Such a fragment is expected to be embedded within 



a collapsing envelope; thus, its mass grows with time. Al- 
though these fragments later undergo another phase of col- 



lapse due to the dissociation of molecular hydrogen (Larson 



196E) in order to become stars (or brown dwarfs), they are 



unlikely to sub-fragment (e.g. Boss 1989; Bate 1998, 2002). 
Thus, the opacity li mit sets a minimum fragment mass of 
« 10 Mj ( |Boss 1988j ). 

To model the opacity limit for fragmentation without 
performing full radiative transfer, we use a gas equation of 
state given by p — Kp v , where K \s a, measure of the entropy 
of the gas. The value of rj varies with density as 



V = 



1, p < 10 -13 g cm -3 , 
7/5, p > 10~ 13 g cm" 3 . 



(1) 



We take the mean molecular weight of the gas to be p = 
2.46. The value of K is defined such that when the gas is 
isothermal K = c%, with c s = 1.84 x 10 4 cm s" 1 at 10 K, 
and the pressure is continuous when the value of r\ changes. 

This equation of state reproduces the temperature- 
density relation of molecular gas during spherically- 
symmetric collapse (as calculated with frequency-dependent 
radiative transfer) to an accuracy of better than 20% i n the 
non-isothermal regime u p to densities of 10 -8 



g cm 3 ([Ma- 



sunaga fc Inutsuka 2000 ). Thus, our equation of state should 



model collapsing regions well, but may not model the equa- 
tion of state in protostellar discs particularly accurately due 
to their departure from spherical symmetry. 



2.2 Sink particles 

The opacity limit results in the formation of distinct 
pressure-supported fragments. As these fragments accrete, 
their central density increases, and it becomes computation- 
ally impractical to follow their internal evolution due to the 
short dynamical time-scales involved. Therefore, when the 
central density of a fragment exceeds 10 -11 g cm -3 , we in- 
sert a sink particle into the calculation (Bate et al. 1995). 
Gas within radius r acc = 5 AU of the fragment's centre (i.e. 
the location of the SPH particle with the highest density) is 
replaced by a point mass with the same mass and momen- 
tum. The fragment is only replaced if it is gravitationally 
bound. Any gas that later falls within this radius is accreted 
by the point mass if it is bound. Thus, we can only resolve 
discs around sink particles if they have radii <; 10 AU. 

Since all sink particles are created from the pressure- 
supported fragments, their initial masses are ~ 10 Mj, as 
given by the opacity limit for fragmentation. Subsequently, 
they may accrete large amounts of material to become 
higher-mass brown dwarfs (^ 75 Mj) or stars (<; 75 Mj). 

The gravitational acceleration between two sink parti- 
cles is Newtonian for r > 4 AU, but is softened within this 



radius using spline softening (Benz 199C). The maximum 
acceleration occurs at a distance of « 1 AU; therefore, this 
is the minimum binary separation that can be resolved. 



2.3 Initial conditions 

The initial conditions consist of a large-scale, turbulent 
molecular cloud. The cloud is spherical and uniform in den- 
sity with a mass of 50 Mq and a diameter of 0.375 pc. At 
the temperature of 10 K, the mean thermal Jeans mass is 
1 Mq (i.e. the cloud contains 50 thermal Jeans masses). An 
initial supersonic turbulent velocity field is imposed on the 
cloud in the same manner as Ostriker, Stone & Gammie 
( fcOO l| ). We generate a divergence-free random Gaussian ve- 
locity field with a power spectrum P(k) tx k~ 4 , where k 



is the wavenumber. In three dimensions, this gives a veloc- 
ity dispersion that varies with distance, A, as cr(A) oc A 1 / 2 
in agreement with the observed Larson scaling relations for 
molecular clouds (Larson 1981). The velocity field is nor- 



malised so that the kinetic energy of the turbulence equals 
the magnitude of the cloud's gravitational potential energy. 



2.4 Resolution 

The local Jeans mass must be resolved throughout the cal- 
culation to correctly model fragmentation (Bate & Burkert 
199 7; Tru elove et al. 1997; Whitworth 1998). Bate & Burk- 
ert (1997) found that this requires <; 2A?neigh SPH particles 
per Jeans mass; -/V 110 i g h is insufficient. We have repeated their 
calculation using different numbers of particles and find that 
1.5iV no i g h = 75 particles is also sufficient to resolve fragmen- 
tation (Bate, Bonnell & Bromm, in preparation). The min- 
imum Jeans mass in the calculation presented here occurs 
at the maximum density during the isothermal phase of the 
collapse, p = 10" 13 g cm" 3 , and is « 0.0011 M© (1.1 Mj). 
Thus, we use 3.5 x 10 6 particles to model the 5O-M0 cloud. 
The calculation required approximately 95000 CPU hours 
on the SGI Origin 3800 of the United Kingdom Astrophys- 
ical Fluids Facility (UKAFF). 
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Figure 1. Top panels: Time sequence showing the formation of two brown dwarfs via disc fragmentation. Spiral density waves in a 
gravitationally-unstable circumbinary disc (panel 1; t = 5160 years) fragment to form a third star (panel 2; t = 5735 years) and two 
brown dwarfs (panel 3; t = 6020 years). The resulting multiple system undergoes chaotic evolution. The first brown dwarf (box; red 
lines) is eventually ejected from the system, while the second is put into a highly eccentric orbit (panel 4; t = 10490 years). Much later 
the second brown dwarf (triangle; blue lines) is also ejected (panel 5; t = 17340 years) leaving a hierarchical triple system. Each panel is 
600 AU across. The colour scale shows the logarithm of column density ranging from 3 x 10" 1 to 3 x 10 3 g cm . Lower panels: The 
trajectories and masses of the five objects during the ejection of the two brown dwarfs. The left panel has the same orientation as the 
time sequence. The centre panel shows the third dimension. The right panel shows the evolution of the masses of the objects. When the 
brown dwarfs are ejected from the dense gas accretion on to them is halted and they are consequently unable to attain stellar masses. 
Times are given in years from the beginning of star formation in the cloud. 



Figure 2. Time sequence showing the formation of two brown dwarfs via the fragmentation of collapsing filaments of molecular gas. 
Two brown dwarfs (box, triangle; red and blue lines) fragment out of filaments (panels 1, 2, & 3; t = 14300, t = 14870 and t = 15440 
years). Both fall towards an existing multiple system, one almost radially along its filament (triangle). The first brown dwarf (square; red 
line) is quickly ejected from the multiple system (panels 3 & 4; t = 15440 and t = 17150 years). Much later, after several crossing times 
of the multiple system, the second brown dwarf (triangle; blue line) is also ejected (panel 5; t = 23710 years). Each panel is 1000 AU 
across. The colour scale shows the logarithm of column density ranging from 3 X 10" 1 to 3 X 10 3 g cm -2 . Lower panels: The trajectories 
and masses of the objects during the ejection of the two brown dwarfs. The left panel has the same orientation as the time sequence. 
The centre panel shows the third dimension. The right panel shows the evolution of the masses of the objects. When the brown dwarfs 
are ejected from the dense gas, accretion onto them is halted and they are consequently unable to attain stellar masses. Times are given 
in years from the beginning of star formation in the cloud. For clarity, only the trajectories of the two brown dwarfs are coloured. 



3 RESULTS 

3.1 Evolution of the cloud 

The hydrodynamical evolution of the cloud produces shocks 
which decrease the turbulent energy that initially supported 
the cloud. In parts of the cloud, gravity begins to dominate 
and dense self-gravitating cores form and collapse. These 
dense cores are the sites where the formation of stars and 
brown dwarfs occurs. In all, 23 stars and 18 brown dwarfs 
are produced. An additional 9 objects have substellar masses 
when we stop the calculation, but are still accreting. Three 
of these are likely to remain substellar, but the other six 
would probably become stars. The cloud's evolution and the 
properties of the stars and brown dwarfs will be discussed 
in detail in subsequent papers. In this letter, we concentrate 
on the mechanism by which the brown dwarfs form. 



3.2 The formation of brown dwarfs 

Although the thermal Jeans mass is initially 1 Mq , the cloud 
produces roughly equal numbers of brown dwarfs and stars. 
By what mechanism(s) do these brown dwarfs form? 

As mentioned in Section 2.2, all the stars and brown 
dwarfs begin as fragments whose initial masses are given by 
the opacity limit for fragmentation. Thus, all the objects be- 
gin with masses of a few Jupiter masses. Those that become 
stars subsequently accrete large quantities of gas from the 
dense cores in which they form, while those that remain as 
brown dwarfs do not. 

Without exception, we find that those objects that end 
up as brown dwarfs begin as opacity-limited fragments in 
dynamically-unstable multiple systems and are ejected be- 
fore they can accrete enough gas to become stars. In no case 
does a dense core collapse to form a single brown dwarf or 
a bina ry brown dwarf system alone. This general mecha- 



sulting properties of the brown dwarfs could only be conjec- 
tured. 

By tracing the evolution of the brown dwarfs in our cal- 
culation back to their fragmentation out of the gas, we find 
that brown dwarfs form in two different ways. Roughly three 
quarters (14 of the 18 definite brown dwarfs) form via the 
fragmentation of gravitationally-unstable circumstellar discs 
(e.g. Figure 1). The discs fragment due to rapid accretion 
(e.g. Bonnell 1994; Bonnell & Bate 1994a,b; Whitworth et 
al. 1995; Burkert, Bate & Bodenheimer 1997) and/or tidal- 
perturbations during stellar encounters (Boffin et al. 1998; 
Watkins et al. 1998a, b). The fragments destined to become 
brown dwarfs then interact with the stars and brown dwarfs 
to which they are bound until they are ejected from the 
dense gas, limiting their masses to be substellar. The re- 
maining brown dwarfs (4 of the 18) originate from the frag- 
mentation of dynamically-collapsing, filamentary molecular 
gas (e.g. Figure 2; Bonnell et al. 1991). The fragments ini- 
tially form on their own, but quickly fall into unstable mul- 
tiple systems and are ejected from the dense gas before they 
can attain stellar masses. This latter mechanism is similar 



to that discussed by Reipurth & Clarke (2001) 



nism h as been proposed by Reipurth fc Clarke ( 2001) and 



Watkins et al. (1998b), but the main processes involved in 
the formation of the unstable multiple systems and the re- 



Although the initial thermal Jeans mass in the cloud is 
1 Mq, the gas in the circumstellar discs and dynamically- 
collapsing gas is at a much higher density. Thus, the local 
Jeans mass and Jeans length are much smaller (both scale 
with density as p _1,/2 ). Consequently, when gravitational 
collapse occurs in such an environment, the opacity-limited 
fragment that forms has a very limited reservoir of gas from 
which to accrete. In order to attain a large mass it must 
accrete most of its mass from distances much greater than 
the Jeans length of the gravitationally unstable clump from 
which it formed. However, since it is in a multiple system it 
must compete for the gas (e.g. Bonnell et al. 1997). Further- 
more, in an unstable multiple system, the mos t frequently 
eject ed object is that having the lowest mass ( van Albada 
196S ) . These two constraints frequently conspire to eject a 
low-mass fragment from the dense gas in which it forms be- 
fore it has accreted enough gas to become a star. 
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If brown dwarfs are produced by ejection from multi- 
ple systems, one might expect the fraction of binary brown 
dwarfs to be low since the binary must be ejected without 
being disrupted. However, a close binary brown dwarf might 
be ejected from a multiple system and avoid disruption if the 
other components of the system were widely separated. 

This calculation does not produce any definite binary 
brown dwarfs, but at the end of the calculation there is a 
close binary brown dwarf within an unstable multiple sys- 
tem. The system consists of 7 objects. The close binary 
brown dwarf (semimajor axis 6 AU) is in a hierarchical bi- 
nary (semimajor axis 84 AU, eccentricity 0.20) with a close 
stellar binary (semimajor axis 7 AU). This quadruple system 
orbits a triple system at 330 AU with an eccentricity of 0.18. 
The triple system consists of a close (9 AU) stellar binary or- 
bited by a brown dwarf companion (semimajor axis 90 AU, 
eccentricity 0.82). This system is dynamically unstable and 
will undergo further evolution. Unfortunately, we are unable 
to follow it any further. Since the binary brown dwarf is very 
close, it is possible that it will survive the dissolution of this 
multiple system. Even so, this would result in only one bi- 
nary brown dwarf system and « 20 single brown dwarfs. We 
therefore conclude that the formation of close binary brown 
dwarfs may be possible, but the fraction of brown dwarfs 
with a brown dwarf companion is very low (^ 5%). 

Observationally, the frequency of brown dwarf binaries 
is not yet clear. Reid et al. (2001) find that 4 out of 20 
brown dwarf primaries have companions giving a fraction 
of ~ 20%. However, this survey is magnitude limited rather 
than volume limited and is therefore likely to overestimate 
the true frequency of brown dwarf binaries. It is interest- 
ing to note that none of the binary brown dwarf systems 
currently known have projected separations > 10 AU (Reid 
et al. 2001). This is consistent with their having survived 
ejection from unstable multiple systems. 

Similarly, the survival of large discs around brown 
dwarfs is not an obvious outcome of the ejection mecha- 
nism of brown dwarf formation. Of the 18 definite brown 
dwarfs, most have close dynamical encounters (^ 20 AU 
for 14 of the 18). Of the remaining four, three are ejected 
before they can accrete gas with the high specific angular 
momentum required to form large discs. Only one brown 
dwarf is ejected with a resolved disc (radius « 60 AU). The 
other brown dwarfs are likely to possess smaller discs, but 
in this calculation we are unable to re solve d iscs with radii 
<, 10 AU (Section 2.2). Muench et al. (|200l|) find » 65% of 



brown dwarfs in the Orion Trapezium cluster have infrared 
excesses indicative of discs. Our results show that if brown 
dwarfs are formed by the ejection mechanism, most of these 
discs should have radii 5j 10 AU. 

Most brown dwarfs in our calculation are produced via 
gravitational instabilities in massive discs. The ease with 
which disc fragmentation occurs depends primarily on the 
rate at which it accretes mass from the surrounding cloud 
(Bonnell 1994) and the disc's equation of state (e.g. Pickett 
et al. 2000). The density of the gas in the discs that frag- 
ment to form brown dwarfs is high enough that the gas is in 
the r\ — 7/5 regime (Section 2.1). Thus, the gas resists frag- 
mentation far more than it would if an isothermal equation 
of state was used, although our equation of state does not 
include heating from shocks. On the other hand, because 
the flattened disc geometry may allow more rapid cooling 



than a spherically-symmetric geometry, real discs may be 
cooler and more unstable than those we model here. Hence, 
the number of brown dwarfs may be even greater. Another 
factor is that stars and brown dwarfs cannot merge in our 
calculation. Mergers may reduce the final number of brown 
dwarfs. In summary, a more definitive prediction will have 
to wait until a large-scale calculation is performed with ra- 
diative transfer and even higher resolution. For the present, 
we have demonstrated that the fragmentation of a turbu- 
lent molecular cloud is capable of forming similar numbers 
of brown dwarfs and stars. This is in agreement with obser- 
vational surveys of brown dwarfs in the solar neighbourhood 
(Reid et al. 1999) and in clusters (e.g. Luhman et al. 2000 
and references therein). 



4 CONCLUSIONS 

We have presented results from the first calculation to 
demonstrate the formation of brown dwarfs from the frag- 
mentation of a turbulent molecular cloud. The calculation 
resolves the opacity limit for fragmentation and follows nu- 
merous brown dwarfs until accretion onto them has ceased. 

We find that the star-formation process produces 
roughly equal numbers of stars and brown dwarfs, in agree- 
ment with recent observations. Examining the mechanisms 
by which the brown dwarfs form, we find that they are the 
result of the ejection of fragments from the dense gas in 
which they form by dynamical interactions in unstable mul- 
tiple systems. This occurs before they are able to accrete to 
stellar masses. Three quarters of the brown dwarfs fragment 
out of gravitationally-unstable discs, while the remainder 
form in collapsing filamentary flows of high-density gas. 

The calculation indicates that close binary brown dwarf 
systems (separations ^ 10 AU) might be able to survive 
the ejection process. However, such systems should be very 
rare (frequency ^ 5%) because of the close dynamical in- 
teractions that are involved in the ejection of brown dwarfs. 
Similarly, large discs (radii <; 10 AU) around young brown 
dwarfs should also be rare (frequency « 5%). Observations 
show that close binary brown dwarfs and circumstellar discs 
surrounding brown dwarfs do exist, but more detailed sur- 
veys are required to test the predictions made here. 



ACKNOWLEDGMENTS 

The computations reported here were performed using the 
UK Astrophysical Fluids Facility (UKAFF). We thank the 
referee for prompting several improvements to the letter. 

REFERENCES 

Basri G., Martin E. L., 1999, AJ, 118, 2460 
Bate M. R., 1998, ApJ, 508, L95 
Bate M. R. , 2002, in preparation 

Bate M. R., Bonnell I. A., Price N. M., 1995, MNRAS, 277, 362 
Bate M. R., Burkert A., 1997, MNRAS, 288, 1060 
Benz W., 1990, in Buchler J. R., ed., The Numerical Model- 
ing of Nonlinear Stellar Pulsations: Problems and Prospects. 
Kluwer, Dordrecht, p. 269 
Benz W., Bowers R. L., Cameron A. G. W., Press W., 1990, ApJ, 
348, 647 



The formation of brown dwarfs 5 

Boffin H. M. J., Watkins S. J., Bhattal A. S., Francis N., Whit- 
worth A. P., 1998, MNRAS, 300, 1189 

Bonnell I. A., 1994, MNRAS, 269, 837 

Bonncll I. A., Bate M. R., 1994a, MNRAS, 269, L45 

Bonnell I. A., Bate M. R., 1994b, MNRAS, 271, 999 

Bonncll I. A., Bate M. R., Clarke C. J., Pringle J. E., 1997, MN- 
RAS, 285, 201 

Bonnell I., Martel H., Bastion P., Arcoragi J. -P., Benz W., 1991, 

ApJ, 377, 553 
Boss A.P., 1988, ApJ, 331, 370 
Boss A.P., 1989, ApJ, 346, 336 

Burkert A., Bate, M.R., Bodenheimcr P., 1997, MNRAS, 289, 497 
Johnstone D., Wilson C. D., Moriarty-Schieven G., Joncas G., 

Smith G., Gregersen, E., Fich, M., 2000, ApJ, 545, 327 
Larson R. B., 1969, MNRAS, 145, 271 
Larson R. B., 1981, MNRAS, 194, 809 

Larson R. B., 1999, in Nakamoto T., ed., Star Formation 1999. 

Nobeyama Radio Observatory, p. 336 
Low C, Lyndcn-Bcll D., 1976, MNRAS, 176, 367 
Luhman K. L., Rickc G. H., Young E. T., Cotcra A. S., Chen H., 

Rieke M. J., Schneider G., Thompson R. I., 2000, ApJ, 540, 

1016 

Masunaga H., Inutsuka S., 2000, ApJ, 531, 350 
Monaghan J. J., Gingold R. A., 1983, J. Comput. Phys., 52, 374 
Motte F., Andre P., Neri R., 1998, A&A, 336, 150 
Muench A. A., Alvcs J., Lada C. J., Lada E. A., 2001, ApJ, 558, 
L51 

Nakajima T., Oppenheimer B. R., Kulkarni S. R., Golimowski D. 

A., Matthews K., Durrance S. T., 1995, Nature, 378, 463 
Ostriker E. C, Stone J. M., Gammic C. F., 2001, ApJ, 546, 980 
Pickett B. K., Duriscn R. H., Casscn P., Mejia A. C, 2000, ApJ, 

540, L95 

Rebolo R., Zapatero-Osorio M. R., Martin E. L., 1995, Nature, 
377, 129 

Rees M. J., 1976, MNRAS, 176, 483 

Reid I. N., Kirkpatrick J. D., Liebert J., Burrows A., Gizis J. E., 

Burgasser A., Dahn C. C, Monet D., Cutri R., Bcichman C. 

A., Skrutskie M., 1999, ApJ, 521, 613 
Reid I. N., Gizis J. E., Kirkpatrick J. D., Koerner D. W., 2001, 

AJ, 121, 489 
Reipurth B., Clarke, C, 2001, AJ, 122, 432 

Truelove J. K., Klein R. I., McKee C. F., Holliman J. H.II, Howell 
L. H., Grccnough J. A., 1997, ApJ, 489, L179 

van Albada T. S., 1968, Bull. Astron. Inst. Ncth., 19, 479 

Watkins S. J., Bhattal A. S., Boffin H. M. J., Francis N, Whit- 
worth A. P., 1998a, MNRAS, 300, 1205 

Watkins S. J., Bhattal A. S., Boffin H. M. J., Francis N, Whit- 
worth A. P., 1998b, MNRAS, 300, 1214 

Whitworth A. P., 1998, MNRAS, 296, 442 

Whitworth A. P., Chapman S. J., Bhattal A. S., Disney M. J., 
Pongracic H., Turner J. A., 1995, MNRAS, 277, 727 



This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0206365vl 



This figure "fig2.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0206365vl 



